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ABSTRACT 

The  specific  conductance  of  three  sea  water  samples  (33.  634°/oo, 
34.099°/oo  and  35.001°/oo)  was  measured  at  twelve  temperatures  from 
0°C  to  21.00°C  and  at  nine  pressures  from  atmospheric  to  4,141.8 
decibars  above  atmospheric.     Polynomials  for  specific  conductance  as 
a  function  of  temperature  and  as  a  function  of  pressure   are   computed 
from  the  data,  but  the  temperature-pressure  effect  is  obscure. 

Two  empirical  equations  that  relate  conductivity  ratio  to  temper- 
ature and  pressure  are  also  developed.    These  equations  better  describe 
the  temperature-pressure  effect  on  conductivity  at  the  pressures 
studied  which  correspond  to  the  upper  4,000  meters  of  the  ocean. 
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I.      INTRODUCTION 

In  recent  years,  measurement  of  the  electrical  conductivity  of  sea 
water  has  become  the  most  common  method  of  determining  the  salinity 
of  sea  water.     As  a  result,  UNESCO  has  defined  salinity  in  terms  of  a 
conductivity  ratio  (International  Oceanographic  Tables,  1966).    Salinity 
has  been  so  defined  because  there  is  good  evidence  that  there  is  better 
correlation  between  the  density  and  salinity  calculated  from  conductivity 
than  between  the  density  and  salinity  calculated  from  chlorinity    (Cox 
et  al.  ,  1962). 

Conductivity  ratio  is  the  ratio  of  the  conductivity  of  an  electrolyte 
of  unknown  concentration  to  the  conductivity  of  a  known  standard  con- 
centration of  the  electrolyte.    In  the  case  of  sea  water,  the  standard 
is  Standard  Sea  Water  with  an  exactly  measured  value  near  35°/oo  (parts 
of  dissolved  salts  per  one  thousand  parts  of  sea  water).    The  ratio  is 
formed  with  both  the  sample  and  standard  at  the  same  temperature  and 
at  a  pressure  of  one  atmosphere.    This  ratio  can  be  represented  as  : 
K(S,T,0) /k(35,T  ,0)  where  the  pressure  variable  is  "sea  pressure"  so 
that  zero  denotes  a  pressure  of  one  atmosphere. 

Commonly ,  conductivity  ratio  of  sea  water  is  measured  using  a 
laboratory  salinometer  that  is  calibrated  several  times  a  day,  using 
standard  sea  water,  when  in  use.    The  measured  ratio  is  then  converted 
to  a  ratio  at  the  temperature  (15°C)  for  which  salinity  is  defined  by 
UNESCO. 
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Conductivity  ratio  is  also  measured  using  an  in  situ  salinometer 
which  measures  the  ratio  K(S  ,T  ,P)/K(S  ,T  ,0)  and  converts  this  through 
instrumentation  to  a  ratio  K(S,T  ,P)/K(35,  15,   0) ,  which  in  turn  is  used 
to  obtain  salinity.    The  measurement  of  conductivity  ratio  in  situ    re- 
quires a  knowledge  of  the  effects  of  temperature  and  pressure  on 
conductivity.    The  in  situ  salinometer  has  been  a  major  development  in 
the  collection  of  oceanographic  data  and  has  been  described  by  Hamon 
(1955),  Prichard  (1959)  and  Siedler  (1963). 

The  effect  of  temperature  and  pressure  on  ionic  conductivity  has 
been  studied  previously.    However ,  the  temperature-conductivity  re- 
lationship has  usually  been  studied  at  atmospheric  pressure ,  and  the 
pressure-conductivity  relationship  at  only  a  few  temperatures. 

With  the  increased  use  of  in  situ  measurements  in  oceanography , 
this  study  was  undertaken  in  an  effort  to  investigate  further  the 
temperature-pressure  relationship  with  conductivity. 

II.      PREVIOUS  INVESTIGATIONS 

Many  studies  on  the  conductivity  of  aqueous  solutions  have  been 
made  since  the  turn  of  the  century ,  but  only  in  the  last  thirty  years 
has  serious  research  using  sea  water  been  carried  out.    The  majority 
of  the  investigations  to  date  have  concentrated  on  temperature  effects , 
and  only  since  the  late  1950's  have  studies  of  pressure  effects  become 
prominent. 


A.     STUDIES  MADE  AT  ATMOSPHERIC  PRESSURE 

Thomas,  Thompson,  and  Utterback  (1934)  conducted  one  of  the 
earliest  and  most  extensive  studies  of  the  effect  of  temperature  on 
the  conductivity  of  sea  water.    They  studied  sea  water  that  ranged  in 
chlorinity  from  l°/oo  to  22°/oo  at  temperatures  of  0°C ,  5°C ,  10°C , 
20°C  and  25°C.      The  results  of  this  study  remained  for  many  years 
the  best  information  available ,  and  was  adequate  for  salinity  deter- 
minations using  inductive  or  conductive  salinometers  where  specific 
conductance  of  a  sample  is  compared  to  a  standard  by  forming  a 
conductivity  ratio. 

Pollack  (1954)  criticized  the  work  of  Thomas  et  al.  because  of 
experimental  errors  from  two  sources:      (1)  Thomas  et  al.  used  the 
results  of  Parker  and  Parker  (1924)  for  the  specific  conductances  of 
standard  KC1  solutions  to  calibrate  the  conductivity  cell.    The  Parker 
and  Parker  values  are  in  error  due  to  capacitive  effects  between  parts 
of  the  cell  that  were  of  opposite  polarity  according  to  Jones  and 
Bradshaw  (1933).      Jones  and  Bradshaw  developed  new  values  for  specific 
conductance  of  standard  KC1  solutions.         (2)  the  conductivity  cells 
(Washburn  pipet  type)  used  by  Thomas  et  al.  allowed  a  capacitive  shunt 
(similar  to  that  of  Parker  and  Parker)  that  could  produce  an  error. 

These  two  sources  of  error  are  opposite  in  sign  and  would  tend  to 
cancel  each  other,  but  to  what  extent  is  indeterminate. 

Weyl  (1964)  reviewed  the  data  of  Thomas  et  al.    and  developed  a 
polynomial  for  evaluating  the  change  in  specific  conductance  ( ^- ) 
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at  atmospheric  pressure  as  a  function  of  chlorinity  and  temperature. 
Weyl  also  developed  a  polynomial  for  specific  conductance  (K)  as  a  function 
of  chlorinity  and  temperature. 

Malmberg  (1965)  has  studied  the  conductivity  of  very  dilute  solutions 
of  sea  water  such  as  would  result  from  a  distillation  process. 

Reeburgh  (1965)  measured  the  specific  conductance  of  Red  Sea  water 
that  had  been  diluted  to  various  chlorinities  at  several  temperatures. 
Reeburgh 's  results  are  in  fair  agreement  with  those  of  Thomas  et  al. 
when  expressed  as  error  in  "equivalent  chlorinity'.' 

Park  (1964a)  showed  that  dissolution  of  CaCC>3  in  sea  water  may  re- 
sult in  as    much  as  0.006°/oo  change  in  salinity  if  conductivity  measure- 
ments are  used  to  estimate  salinity.     Park  (1964b)  also  determined  the 
partial  equivalent  conductances  of  16  of  the  electrolytes  in  sea  water. 

B.     STUDIES  MADE  AT  ELEVATED  PRESSURES 

Adams  and  Hall  (1931)  conducted  one  of  the  earliest  studies  of  the 
effect  of  pressure  on  conductivity.     They  made  measurements  on  various 
concentrations  of  NaCl  solutions  and  on  a  few  other  prepared  solutions. 

Hamon  (1958),  apparently,  was  the  first  to  investigate  the  effect 
of  pressure  on  sea  water.    The  study  occurred  during  his  development 
of  an  in  situ  salinometer.    Measurements  were  made  on  a  single  sample 
of  artificial  sea  water  at  four  different  temperatures.    Pressure  was 
increased  in  three  steps  by  using  compressed  nitrogen  up  to  a  pressure 
of  about  900  decibars.     Adiabatic  compression  was  assumed  and  an 
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appropriate  temperature  correction  applied  after  each  pressure  increase. 
Hamon  described  his  results  as  only  tentative  with  an  estimated  accuracy 
of  the  pressure  coefficients  of  5%.     He  concluded  that  the  relationship 
between  pressure  and  conductivity  was  linear. 

Home  and  Frysinger  (1963)  conducted  an  investigation  of  the  proper- 
ties of  sea  water  at  elevated  pressures,  including  electrical  conductivity. 
They  measured  the  specific  conductance  of  three  samples  of  artificial 
sea  water  at  0°C ,  5°C,  15°C  and  25°C  at  pressures  up  to  13,800  decibars 
(approximately  a  depth  of  13,800  meters)  in  increments  of  1,380  deci- 
bars.   Temperature  was  controlled  by  a  constant  temperature  bath  to 
+0.  05°C.    Home  and  Frysinger  found  that  conductivity  as  a  function  of 
pressure  was  linear  up  to  6,890  decibars  and  that  the  slope  of  the 
conductivity  vs.  pressure  plot  was  independent  of  temperature  in  the 
range  studied.    They  developed  a  linear  equation  that  relates  the  specific 
conductance  at  any  pressure ,  temperature  and  chlorinity  to  the  specific 
conductance  at  atmospheric  pressure ,  and  at  the  same  temperature  and 
chlorinity. 

Bradshaw  and  Schleicher  (1965)  made  measurements  of  the  fractional 
increase  in  electrical  conductance  on  three  different  solutions  of 
salinity:    31°/oo,  33°/oo  and  39°/oo.    Their  measurements,  at  six 
temperatures  and  pressures  up  to  10,338  decibars,  were  the  first  to 
be  made  using  natural  sea  water.    The  pressure  was  increased  in  steps 
of  1,723  decibars  for  each  sample  at  each  temperature.    They  estimated 
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that  the  measurements  had  a  precision  of  better  than  0.  01°/oo    in 
equivalent  salinity  change.      Bradshaw  and  Schleicher  developed  an  em- 
pirical relationship  for  the  per  cent  increase  in  conductance  that  fit 
their  data  with  a  standard  deviation  of  0.002°/oo  in  equivalent  salinity 
error.    However ,  there  is  a  significant  residual  error  near  atmospheric 
pressure. 

Mays  (1968)  made  the  first  study  of  the  effect  of  pressure  in 
which  the  pressure  was  increased  in  relatively  small  increments.    Mays 
studied  five  solutions  (two  sea  water)  at  four  temperatures  up  to  a 
pressure  of  2,000  decibars.     Measurements  were  made  in  a  pressure 
vessel  surrounded  by  a  cooling  coil ,  but  no  provision  was  made  for 
measuring  temperature  of  the  water  samples  above  atmospheric  pres- 
sure.   To  allow  for  heat  generated  by  the  nearly  adiabatic  compression , 
Mays  waited  an  arbitrary  15  minutes  after  each  pressure  increase  before 
recording  any  data.    Mays  concluded  that  the  specific  conductance  is 
related  by  a  higher  order  relationship  to  pressure.    Due  to  the  lack  of 
temperature  measurements,  the  nonlinearities. suggested  by  Mays  may 
be  due  to  temperature  fluctuations. 

Ettle  (1969)  investigated  the  effect  of  pressure  on  conductance  of 
two  natural  sea  water  samples  at  six  temperatures  and  pressures  to 
4,000  decibars.    Ettle  also  estimated  the  effects  that  cause  the  in- 
crease in  specific  conductance  of  sea  water  at  elevated  pressures.    He 
found  that  the  change  in  specific  conductance  with  pressure,  while 
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nearly  linear,  must  be  described  by  a  higher  order  polynomial.    The 
standard  deviation  for  second  and  third  degree  equations  was  0.  041°/oo 
equivalent  error  in  salinity  and  fit  the  Bradshaw  and  Schleicher  equation 
with  a  standard  deviation  of  0.  066°/oo.    Ettle  thus  confirmed  that  the 
higher  order  results  of  Bradshaw  and  Schleicher  are  more  descriptive 
of  the  relationship  between  specific  conductance  and  pressure.    How- 
ever, more  investigation  into  the  effect  of  pressure  at  the  temperatures 
characteristic  of  the  ocean  water  column  is  needed  in  order  to  improve 
the  data  obtained  from  in  situ  measurements  of  salinity. 

III.      EQUIPMENT 

The  equipment  used  in  this  study  was  the  same  as  that  used  by 
Ettle  (1969).    However,  the  pressure  vessel  was  remachined,  the 
conductivity  cell  was  rebuilt,  another  pressure  gauge  was  used  and  an 
oscilloscope  was  utilized  vice  earphones  as  a  bridge  balance  detector. 

A.      PRESSURE  VESSEL 

The  pressure  vessel  was  made  in  three  major  pieces:    the  base, 
the  barrel  and  the  piston    (Figure  1).      The  interior,  where  the  con- 
ductivity cell  was  located ,  of  the  vessel  when  assembled  was  about 
1.  5  inches  in  diameter  and  7  inches  deep. 

The  base  had  a  1.  5  inch  hole  bored  to  a  depth  of  one  inch,  to  serve 
as  a  receptacle  for  the  conductivity  cell  that  extended  below  the  level 
of  the  cell  electrodes.    There  were  three  radial  grooves  cut  into  the 
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Figure  1.       Schemati c  Drawi ng  of  the   Conductivity  Cell    and  Pressure 
Vessel.       A — Conductivity  Cell  ,  B — Glass  Cover,  C — Pressure  Vessel 
Barrel  ,  D — Pressure  Vessel  Piston,  E — Pressure  Vessel  Base. 
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top  of  the  base  two  grooves  which  were  diametrically  opposite  each  other 
and  provided  for  the  cell  external  connecting  leads.     The  third  groove 
provided  the  location  of  the  copper-constantan  thermocouple.    All  leads 
were  wrapped  with  TEFLON  tape  to  provide  electrical  insulation  and 
make  a  pressure  seal. 

The  barrel  of  the  vessel  was  about  6  inches  high  with  a  1.  5  inch 
bore  and  5  inch  outside  diameter.    A  circular  groove  was  cut  into  the 
bottom  of  the  barrel  to  accept  a  neoprene  0-ring  that  made  a  pressure 
seal  when  the  barrel  was  bolted  to  the  base.    The  base  and  barrel  were 
bolted  together  by  six  stainless  steel  bolts. 

The  piston  was  made  slightly  smaller  in  diameter  than  the  bore  of 
the  barrel.  A  pressure  seal  was  obtained  by  using  a  neoprene  O-ring  in 
a  groove  located  near  the  bottom  of  the  piston. 

After  the  vessel  was  assembled,  short  lengths  of  plastic  "spaghetti" 
were  placed  over  each  lead  that  penetrated  the  pressure  vessel ,  then 
two-inch  lengths  of  stainless  steel  tubing  were  installed  over  the 
"spaghetti"  and  secured  to  the  side  of  the  pressure  vessel.    This 
arrangement  provided  mechancical  strength  and  isolated  the  leads 
from  the  cooling  coil. 

The  cooling  coil,  0.  5  inch  copper  tubing,  was  placed  around  the 
assembled  base  and  barrel. 

B.     CONDUCTIVITY  CELL 

The  conductivity  cell  was  constructed  from  a  glass  cylinder  closed 
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at  one  end  and  fitted  with  a  glass  cover  at  the  other  end  (Figure  1).     Two 
glass  stand-offs  were  located  near  the  top  of  the  cell.    These  stand- 
offs held  the  cell  in  position  inside  the  pressure  vessel.    Platinum  foil 
electrodes  approximately  one  square  centimeter  in  area  were  located 
about  one  inch  from  the  bottom  of  the  glass  cylinder  and  tack-welded 
to  platinum  wire  leads  that  penetrated  the  cylinder  wall.     The  wire 
leads  provided  the  necessary  external  connections  for  measuring  con- 
ductance.   The  platinum  leads  were  fixed  to  the  cylinder  using  PLIOBOND 
electrical  sealant. 

The  foil  electrodes  were  coated  with  platinum  black  using  Braun 
Chemical  Co.  platinizing  solution  and  a  3  volt  battery.    After  being 
platinized,  the  cell  was  kept  full  of  distilled  water  when  not  in  use. 

C.  PRESSURE 

The  assembled  pressure  vessel  with  piston  in  place ,  was  put  in  a 
hydraulic  press  and  pressure  applied  to  the  base  and  top  of  the  piston. 
Hydraulic  pressure  was  measured  using  an  eighteen  inch  Heise  gauge 
(H3156)  with  a  range  to  6/300  psi.      The  pressure  inside  the  vessel  could 
then  be  computed  since  the  measured  hydraulic  pressure  is  directly 
related — by  a  ratio  dependent  upon  the  area  of  the  respective  pistons — 
to  the  pressure  inside  the  vessel. 

D.  TEMPERATURE 

Temperature  in  the  vessel  was  controlled  by  pumping  an  ethylene 
glycol-water  mixture  from  a  Formatemp  (2075B)  thermostatically 
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controlled  bath  through  the  cooling  coil  that  surrounded  the  vessel.     The 
bath  could  control  temperature  to  within  +  0.02°C.      The  temperature 
inside  the  vessel  was  measured  with  a  copper-constantan  thermocouple 
located  at  the  same  level  as  the  platinum  foil  electrodes.    Potential 
was  measured  by  a  Leeds  and  Northrup  precision  millivolt  potentiometer 
(8690). 

E.  CONDUCTIVITY 

The  resistivity — the  inverse  of  conductivity — of  the  water  samples 
in  the  cell  was  measured  with  a  Wheatstone  bridge.    The  internal  decade 
resistors  could  be  adjusted  in  one  ohm  steps.    In  order  to  obtain  the 
accuracy  required,  three  decade  resistance  boxes  (adjustable  in  0. 1  ohm 
steps)  were  connected  in  parallel  and  connected  to  the  bridge.    Bridge 
balance  was  observed  using  a  Heathkit  oscilloscope  (EUW-25). 

F.  EQUIPMENT  ACCURACY  AND  PRECISION 

Pressure  could  be  measured  to  an  accuracy  of  5  psig  and  the  piston 
arrangement  increased  the  pressure  in  the  vessel  by  a  factor  of  about 
1.  5  thus  the  pressure  inside  the  vessel  was  accurate  to  7.  5  psi.    This 
corresponds  to  an  accuracy  of  5.  2  decibars  or  an  equivalent  salinity 
error  of  +0.  003  °/oo. 

Temperature  could  be  controlled  to  +0.  02°C  and  measured  with  an 
accuracy  of  +0.01°C.      This  accuracy  in  temperature  corresponds  to  an 
equivalent  salinity  error  of  +0.02°/oo. 
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Resistance  of  the  cell  could  be  measured — using  the  external  decade 
resistors — with  an  accuracy  of  0.01  ohm.      Since  resistances  measured 
were  of  order  25  ohms,  and  the  cell  constant  was  about  0.7,  this 
accuracy  in  resistance  is  equivalent  to  an  accuracy  of  +0.  00001  mho/cm 
in  specific  conductance  or  an  equivalent  salinity  error  of  +0.01°/oo. 

Experimental  accuracy  and  precision  are  summarized  in  the  following 
table : 


Variable 


Accuracy       Precision 


Equivalent 
Salinity 
Error  (°/oo) 


Pressure  (decibars) 

Temperature  (°C) 

Specific  Conductance 
(mho /cm) 


5.2 

5.2 

+0.003 

0.01 

0.02 

+0.02 

0.00001 

0.00001 

+0.01 

Thus ,  the  maximum  total  experimental  error  is  not  greater  than 
+0.033°/oo  error  in  equivalent  salinity,  and  the  most  probable  error 
is  +0.  016°/oo. 

IV.      EXPERIMENTAL  PROCEDURE 

The  conductance  cell  was  calibrated  while  installed  in  the  pressure 
vessel  at  atmospheric  pressure  and  0°C  with  0. 1M  standard  KC1  solution. 
Conductivity  of  three  water  samples — two  natural  sea  water  and 
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standard  sea  water — was  then  determined  at  twelve  temperatures  at 
pressure  intervals  of  500  psi  to  a  maximum  of  4,000  psi. 

A.     CELL  CALIBRATION 

To  accurately  measure ,  from  first  principles ,  the  specific  con- 
ductance of  an  electrolyte ,  such  as  sea  water ,  requires  a  cell  with 
electrodes  of  accurately  determined  area  and  spacing.    However ,  an 
effective  value  of  spacing/area  can  be  determined  by  using  standard 
solutions  of  known  concentration  and  known  specific  conductance.    This 
factor — spacing/area  (L/A) — is  different  for  each  cell  and  is  called 
the  cell  constant.    A  cell  constant  for  any  given  conductivity  cell  can 
be  determined  by: 

K0C  =  KT  (1) 

Where: 

K0     =  observed  conductivity  of  standard  solution  as 
measured  by  a  conductivity  bridge  (mho) 

C       =  L/A  =  cell  constant  (cm     ) 

KT     =  specific  conductance  of  standard  solution  of 
known  concentration  at  a  given  temperature 
from  literature  (mho/cm) 

Once  a  cell  constant  has  been  determined,  the  specific  conductance 

of  an  unknown  electrolyte  can  then  be  calculated  from: 

Kt=  K0C  (2) 
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Where:  K»    =  specific  conductance  of  sample  at  a  known 

temperature  and  atmospheric  pressure 
(mho /cm) 

Kc  =  observed  conductivity  of  sample  as  measured 
by  a  conductivity  bridge  (mho) 

C    =  cell  constant  (cm     ) 

In  this  study,  since  the  platinum  leads  of  the  cell  were  held  in  place 

by  the  pressure  vessel  and  the  cell  penetrations  were  flexible ,  no  changes 

in  the  cell  constant  with  pressure  were  expected.    However,  a  new  cell 

constant  determination  was  necessary  whenever  the  pressure  vessel  and 

cell  were  disassembled. 

B.      STANDARD  SOLUTIONS 

Two  different  standard  solutions  of  0.1M  KC1  were  made  using 
reagent  grade  KC1  crystals  using  the  values  and  procedures  according 
to  Kortum  (1965).    The  KC1  crystals  were  first  dried  by  heating  for 
several  hours  and  allowed  to  cool  in  a    desiccator  before  weighing. 
Duplicate  measurements  of  conductivity  were  made  on  each  solution. 
The  values  of  the  measurements  were  within  0.00002  mho/cm  of  each 
other  and  the  value  used  in  determining  the  cell  constant  was  the  aver- 
age of  these  four  measurements.    The  distilled  water  used  in  making 
the  standard  solutions  was  found  to  have  a  specific  conductance  of 
0.00016  mho/cm  and  was  corrected  for  when  calculating  the  cell 
constant. 
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C.  WATER  SAMPLES 

A  large  parcel  of  natural  sea  water  was  collected  from  the  southern 
end  of  Monterey  Bay.     This  sample  was  divided  into  two  parts;  the  first 
was  used  as  collected;  the  other  was  left  exposed  to  the  air  and  allowed 
to  evaporate  at  room  temperature  until  a  higher  salinity  was  obtained. 
The  salinity  of  each  sample  was  determined  using  a  Hytech  (622)  inductive 
salinometer  and  the  final  salinities  were  33.  634°/oo  and  34.  099°/oo. 
The  third  sample  used  was  several  ampoules  of  standard  sea  water  with 
a  chlorinity  of  19.  3745° /oo. 

For  each  water  sample ,  measurements  were  taken  at  twelve 
temperatures  from  near  0  C  to  21°C  with  most  measurements  below 
10°C.    Measurements  were  taken  in  500  psi  increments  up  to  4,000  psi. 

D.  MEASUREMENTS 

For  each  water  sample  studied,  both  the  conductivity  cell  and 
pressure  vessel  were  filled.    Before  filling,  the  equipment  was  rinsed 
twice  with  distilled  water  and  twice  with  the  particular  water  sample 
to  be  used.    No  aliquot  of  a  sample  was  ever  reused. 

After  filling  the  sample  chamber ,  the  glass  cover  and  piston  were 
installed    and  the  pressure  vessel  placed  in  the  hydraulic  press.     Con- 
ductivity of  the  sample  was  then  determined  at  pressure  increments 
of  500  psi  from  atmospheric  to  4,000  psi  and  from  4,000  psi  back  to 
atmospheric.    The  measurements  at  each  pressure  were  then  averaged 
since  the  readings  at  each  pressure  were  within  0.02  ohm  of  each  other. 
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The  temperature  bath  was  adjusted  to  give  a  desired  temperature 
within  the  vessel,  and  after  each  pressure  change  the  bridge  was  balanced 
and  temperature  measured  alternately  until  three  consistent  readings 
were  obtained  (within  0.01  ohm  resistance  and  0.02°C);  then  pressure, 
temperature  and  bridge  resistance  were  recorded.    Since  no  specific 
temperature  was  required,  no  adjustment  of  the  bath  for  adiabatic 
heating  or  cooling  was  necessary  following  each  pressure  change. 

All  measurements  were  made  at  a  frequency  of  2,000  Hz,  and 
occasionally  a  value  was  checked  at  1,000  Hz  to  insure  observations 
were  independent  of  frequency. 

V.      EXPERIMENTAL  RESULTS  AND  ANALYSIS 

Experimental  measurements  were  made  on  three  sea  water  samples 
at  twelve  temperatures.    At  each  temperature,  three  sets  of  data  for 
values  of  conductivity  at  nine  pressures  were  obtained.    The  results  are 
tabulated  in  Tables  I  through  III  of  Appendix  A. 

A.      CONVERSIONS 

Before  detailed  analysis  of  the  experimental  data  could  be  performed, 
several  conversions  and  determinations  were  required.    The  raw  data 
consisted  of  gauge  pressures,  sample  resistances,  and  temperatures. 

1.      Pressure 

The  hydraulic  pressure  in  psig  was  converted  to  equivalent 
pressure  inside  the  vessel  and  expressed  in  decibars.      This  estabHshed 
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some  perspective  of  pressure  and  depth  in  the  water  column  since 
1  decibar  corresponds,  approximately,  to  a  depth  of  1  meter. 

The  pressure  inside  the  pressure  vessel  is  related  to  the 
hydraulic  pressure  by  a  ratio  of  the  diameters  of  the  hydraulic  and 
vessel  pistons. 
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PV=Ph\Dv/XF  (3) 

Where:  Py  =  pressure  inside  pressure  vessel  (decibars) 

P^  =  hydraulic  pressure  observed  (psig) 

D^  =  diameter  of  hydraulic  piston  (inches) 

Dy  =  diameter  of  pressure  vessel  piston  (inches) 

F     =  factor  for  converting  psig  to  decibars  =0.  68947 
(Wilson  and  Bradley,  1968) 

2.      Conductivity 

The  conductivity  bridge  measures  the  resistance  of  the  solution 
across  the  electrodes ,  which  can  be  converted  to  specific  conductivity. 

The  resistance  of  the  external  decade  boxes  connected  in 
parallel  was  calculated  by: 


R    =        Rl         x       R2         x        ^3  (4) 

(Rl  X  R2)  +  (R2  X  R3)  +  (Rj_  X  R3) 


Where:     R  =  total  resistance  of  external  decade  boxes  (ohm) 
R-,  =  resistance  of  first  decade  box  (ohm) 
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R2  -  resistance  of  second  decade  box  (ohm) 
R3  =  resistance  of  third  decade  box  (ohm) 

The  total  cell  resistance  was  obtained  by  adding  the  resistance 
measured  by  the  conductivity  bridge,  the  total  external  resistance ,  and 
a  calibration  constant  of  0.06  ohm. 

Rt  =  Rb  +  Re  +  Rc  (5) 

Where:     R^.  =  total  cell  resistance  (ohm) 

R^  =  resistance  measured  by  conductivity  bridge  (ohm) 
Rg  =  total  external  resistance  (ohm) 
Rc  =  calibration  constant  (ohm) 

The  observed  conductivity  in  mhos  is  the  reciprocal  of  the  total 
resistance  (Rt): 


K„  = 


=  _L_  (6) 

Rt 


This  may  then  be  substituted  into  the  expression  for  the  specific  con- 
ductance of  an  unknown  electrolyte ,  to  obtain  the  specific  conductance 
of  an  unknown  electrolyte ,  to  obtain  the  specific  conductance  from  total 
resistance: 


Rt 
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Where:  K$       -  specific  conductance  of  unknown  sample  (mho/cm) 

C        =cell  constant  (cm"   ) 


R.       =  total  measured  resistance  (ohm) 

3.  Temperature 

The  thermocouple  was  calibrated  against  a  precision  mercury 
thermometer  and  sufficient  time  was  permitted  for  thermal  equili- 
bration.   No  conversions  or  corrections  to  temperature  were  required. 
Temperatures  inside  the  vessel,  as  sensed  by  the  thermocouple,  were 
obtained  graphically  from  a  determined  calibration  curve  of  temper- 
ature vs.  potential. 

4.  Salinity 

When  comparing  results  of  this  investigation  with  the  results 
of  previous  studies ,  and  it  was  necessary  to  convert  salinity  to  chlor- 
inity,  the  following  was  used: 

S°/oo  =  1.  80655Cl°/oo  (8) 

The  data  for  all  three  solutions  at  4.  39°C  and  0.  42°C  is  of 
questionable  accuracy  since  this  data  was  taken  first  and  it  was  noted 
that  the  pressure  vessel  piston  had  a  tendency  to  drag  when  pressure 
was  released.    Therefore ,  the  values  for  specific  conductance  at  these 
two  temperatures  may  not  be  at  the  recorded  pressure. 
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B.     DATA  ANALYSIS 

Analysis  of  the  experimental  data  was  performed  using  the  IBM  360 
digital  computer  at  the  Naval  Postgraduate  School.    The  raw  data  was 
punched  onto  IBM  cards  for  processing.    The  computer  program  used 
for  the  analysis  can  be  found  at  the  end  of  this  paper. 

1.      Temperature 

The  coefficients  of  the  least  square  polynomial  up  to  third 
degree ,  for  specific  conductance  as  a  function  of  temperature ,  were 
determined  at  each  pressure  and  are  shown  in  Tables  IV  through  VI 
(Appendix  A).     This  polynomial  is  of  the  form: 


K(S  ,T  ,P)  =  A±  +  A2T  +  A3T2  +  A4T3  (9) 


Where:  A^    =  coefficient  at  each  pressure 

T     =  temperature  in  degrees  centigrade 

The  second  and  third  degree  curves  fit  the  data  with  a  standard 
deviation  of  less  than  0.03°/oo  equivalent  salinity  error  in  all  cases. 
Representative  plots  of  the  second  degree  fit  are  shown  in  Figures  2,3, 
and  4. 

2.      Pressure 

The  coefficients  of  the  least  square  polynomial  up  to  third 
degree ,  were  also  computed  for  specific  conductance  as  a  function  of 
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Figure  2.      Specific  Conductance  vs.  Temperature, 

Salini  ty=33.  634  /oo,  Second  Degree  Curve, 


28 


o 
O 

X 

X 


u 

z 
< 

u 

o 

z 
O 
u 


LU 

CL 


o- 
b 


o- 

b 


1  =4141.8   DECIBARS 

2  ■  1035.5  DECIBARS 


O" 

b 


O" 

b 


O" 

b 


O" 

o 


CI 

O" 

o 


TEMPERATURE     (°C> 


O  0  . 


Figure  3. 


12 


16 


20 


24 


Specific  Conductance  vs.   Temperature, 
Salinity=35.  001°/oo,  Second  Degree  Curve. 
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Specific  Conductance  vs.  Temperature, 
Salinity=34.  099°/oo ,  Second  Degree  Curve. 
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pressure  and  are  shown  in  Tables  VII  through  IX.     This  polynomial  is  of 
the  form: 


K(S,T,P)  =B1  +B2P  +B3P2  +B4P3  (10) 


Where:  B.  =  coefficients  at  each  temperature 

P    =  pressure  in  decibars. 

The  second  and  third  degree  curves  in  this  case  fit  the  data 
with  a  standard  deviation  of  less  than  0.01°/oo  equivalent  salinity  error. 
Representative  plots  of  the  second  degree  fit  are  shown  in  Figures  5,6, 
and  7. 

C.     COMPARISON  WITH  PREVIOUS  INVESTIGATIONS 

The  results  of  this  study  were  compared  with  the  values  of  electrical 
conductivity  of  sea  water  reported  in  previous  investigations.    The 
results  were  compared  in  a  manner  consistent  with  each  previous  study. 
For  example ,  if  the  previous  results  were  reported  as  a  conductivity 
ratio,  the  data  from  this  study  was  converted  to  a  conductivity  ratio 
before  comparisons  were  made. 

la      Atmospheric  Pressure 

The  values  of  specific  conductance  observed  at  atmospheric 
pressure  were  found  to  be  about  1%  less  than  those  reported  by  Thomas 
et  al.   (1934)  or  Reeburgh  (1965).    However,  the  observed  values  are  in 
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Specific  Conductance  vs.  Pressure, 
Salinity=33.  634°/oo ,  Second  Degree  Curve. 
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Figure  6.      Specific  Conductance  vs.   Pressure, 

Salinity=35.  001°/oo,  Second  Degree  Curve. 


T-21.0  0°C 
D 


42  00 


33 


o 

O"  0 


760  1400 

PRESSURE    (oecibars) 


2800 


Figure  7.      Specific  Conductance  vs.  Pressure, 

Salinity=34.  099°/oo ,  Second  Degree  Curve. 
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good  agreement  with  comparable  values  reported  by  Ettle  (1969)  whose 
results  were  also  0.  5%  to  1.  0%  below  those  of  Thomas  et  aL  or 
Reeburgh. 

2.      Elevated  Pressures 

a.  Hamon 

Hamon  (1958)  reported  linear  pressure  coefficients  of  the 
conductivity  ratio  for  sea  water  of  salinity  35.  6°/oo  at  four  temper- 
atures.   The  pressure  range  considered  by  Hamon  was  0  to  1,000  decibars 
above  atmospheric. 

Hamon  made  his  measurements  under  nearly  adiabatic 
conditions  and  then  applied  a  correction  for  the  conductivity  increase 
due  to  adiabatic  heating.    Bradshaw  and  Schleicher  (1965)  suggested  that 
significant  heat  was  probably  lost  and  this  may  explain  why  Hamon' s 
results  are  about  15%  lower  than  the  results  of  this  investigation. 

Further  comparison  with  Hamon's  data  is  not  possible 
since  his  data  shows  increasing  resistance  for  both  increasing  tempera- 
ture and  pressure.    The  reverse  is,  of  course,  true  in  both  cases. 

b.  Home  and  Fry  singer 

Home  and  Frysinger  (1963)  in  their  study  concluded  that  the 
effect  of  pressure  on  conductivity  of  sea  water  is  Hnear  up  to  6,890 
decibars  and  independent  of  temperature  in  the  range  0°C  to  25°C.    How- 
ever ,  as  has  been  previously  shown ,  the  pressure  effect  cannot  accurately 
be  described  as  linear  in  the  range  of  pressures  considered  by  Home  and 
Frysinger. 
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For  purposes  of  comparison,  values  of  specific  conductance 
were  computed  using  the  salinities  and  pressures  of  this  study  in  Home 
and  Frysinger's  equation: 

KP,T,C1=   K1,T   CI  +C  (0.35±0.25)    +    (0.185 ±0.028)C1]10~6P      (11) 

where  pressure  is  in  bars  above  atmospheric  and  chlorinity  (CI)  is  in  parts 
per  mille. 

All  mean  values  of  specific  conductances  computed  using  the  Home 
and  Frysinger  equation  were  within  0.001  mho/cm  of  the  observed  value, 
and  all  observed  values  fall  within  the  rather  large  deviation  described 
by  Home  and  Frysinger. 

c.      Bradshaw  and  Schleicher 

Bradshaw  and  Schleicher  (1965)  essentially  duplicated  the 
work  of  Home  and  Frysinger.    However,  they  reported  their  results  as 
a  polynomial  for  the  fractional  increase  in  conductivity  ratio  of  specific 
conductance  at  pressure ,  P ,  to  the  specific  conductance  at  atmospheric 
pressure.     Their  polynomial  is  first  degree  in  salinity  and  third  degree 
in  temperature  and  pressure.    Their  equation  was  solved  for  the  con- 
ductivity ratio.    The  computed  ratio  was  about  0.  5%  below  the  observed 
conductivity  ratio  calculated  from  the  present  data.     This  corresponds 
to  an  equivalent  salinity  error  of  up  to  0.19°/oo.    This  is  a  significant 
difference ,  but  Anthony  Amos  of  Lamont  Dougherty  Geological  Observa- 
tory observed  at  a  recent  STD  users  conference  sponsored  by  Bissett- 
Berman  Corporation ,  that  in_  situ    salinity  sensors  may  be  in  error  by 
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as  much  as  -0.02  to  0.15°/oo  due  to  apparent  pressure  effects  on  the 
conductivity  head  (Bissett-Berman  Corp.  ,  1971). 

The  conductivity  ratios  calculated  from  the  present  data 
and  the  ratios  computed  using  the  Bradshaw  and  Schleicher  equation 
have  been  plotted  and  representative  plots  are  shown  in  Figures  8,9, 
and  10. 

The  fractional  increase  in  conductivity  ratio ,  when  computed 
using  the  Bradshaw  and  Schleicher  equation ,  should  be  zero  at  atmospheric 
pressure.    However,  examination  of  their  equation  reveals  that  a  residual 
of  at  least  4  x  10      remains.    Thus ,  their  equation  does  not  behave  well 
near  atmospheric  pressure.    This  residual  at  atmospheric  pressure  cor- 
responds to  an  equivalent  salinity  error  of  about  0.  l°/oo. 

In  order  to  avoid  the  problem  of  Bradshaw  and  Schleicher 
of  significant  error  at  pressures  near  atmospheric ,  an  equation  was 
developed  that  more  nearly  fits  the  physical  constraints  on  conductivity 
ratio  at  atmospheric  pressure.    The  equation  developed  has  the  form  of: 

T"1  x  Ln(R)  =  A1  +  A2(P-P0)  +  A2(P-P0)  (12) 

Where:  Aj  =  coefficients  at  each  temperature 

T  =  absolute  temperature  (°K) 

(P-P0)  =  change  in  absolute  pressure  (decibars) 

Ln(R)  =  natural  logarithm  of  the  conductivity 

ratio  K(S,T,P)/K(S,T,P0) 

PQ  =  reference  pressure  of  one  atmosphere 

(10.13  decibars) 
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Figure  8. 


Conductivity  Ratio  vs.   Pressure,  Salinity=33.  634°/oo. 
Points:    Present  data ,  curve:    Bradshaw  and  Schleicher. 
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Figure  9. 


Conductivity  Ratio  vs.  Pressure,  Salinity=35.  001°/oo. 
Points:    Present  data ,  curve:    Bradshaw  and  Schleicher. 
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Figure  10. 


Conductivity  Ratio  vs.  Pressure,  Salinity=34.  099°/oo. 
Points:    Present  data,  curve:    Bradshaw  and  Schleicher. 
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This  equation  still  has  an  ambiguity  at  atmospheric  pressure 
of  e  *        This  is  obviously  a  significant  improvement  if  the  equation 

can  be  verified  by  experiment.     This  equation  fits  the  conductivity 
ratios  calculated  from  the  present  data  with  a  standard  deviation  of 
about  1.  4  x  10"      in  conductivity  ratio  which  corresponds  to  an  equivalent 
salinity  error  of  less  than    +0.01°/oo.     Values  of  the  coefficients  are 
shown  in  Table  X,  XI  and  XII  (Appendix  A). 

Another ,  independent ,  equation  that  fits  the  present  data 
equally  well  was  also  developed.    This  equation  has  the  form: 

(P-P0)~  x  LN  (R)  =B1+  B2T  +  B3T2  (13) 


Where:  B.  =   coefficients  at  each  (P-P0) 

(P-P  )  =    change  in  absolute  pressure  (decibars) 

LN(R)  =   natural  logarithm  of  the  conductivity  ratio 

K(S,T,P)/K(S,T,P0) 

T  =   absolute  temperature  (°K) 


PQ  =   reference  pressure  of  one  atmosphere 

(10.13  decibars) 


Values  of  the  coefficients  for  this  equation  are  shown  in 
Table  XIII,  XIV  and  XV  (Appendix  A). 
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d.  Mays 

Significant  nonlinearities  were  found  by  Mays  (1968)  in  the 
relationship  between  specific  conductance  and  pressure.    He  was  unable 
to  conclude  which  degree  polynomial  to  use  to  describe  the  pressure  effect 
and  none  was  reported.     Consequently,  no  quantitative  comparison  was 
possible. 

Qualitatively,  Mays'  data  showed  large  deviations  from 
linearity  and  occasionally  maxima  and  minima  were  observed.    These 
were  not  observed  in  this  study  nor  have  they  been  reported  in  earlier 
investigations . 

e.  Ettle 

Ettle  (1969)  in  his  investigation  of  the  pressure  effect  on 
conductivity  of  sea  water  reported  a  third  degree  polynomial  in  pressure 
for  conductivity  ratio  that  fit  his  data  with  a  standard  deviation  of 
0.05°/oo  equivalent  salinity  error.    His  data  was  in  fair  agreement  with 
the  results  of  Bradshaw  and  Schleicher. 

Conductivity  ratios  at  three  temperatures  were  computed 
(using  the  polynomial  developed  by  Ettle)  near  corresponding  temperatures 
and  pressures  of  the  present  study.    These  computed  ratios  were  then 
compared  to  conductivity  ratios  calculated  from  the  observed  specific 
conductances.    The  computed  ratios  were  ,  in  all  cases  ,  0.  2%  to  0.  8% 
below  the  observed  ratios.    This  is  not  an  unexpected  result  since  Ettle 's 
equation  has  no  temperature  or  salinity  dependence ,  and  neither  the 
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temperature  nor  salinity  of  the  present  data  was  exactly  the  same  as 
those  used  by  Ettle  when  he  developed  his  equation. 

A  few  randomly  selected  values  of  pressure ,  temperature 
and  specific  conductance  from  Ettle's  data  were  used  in  both  equations 
(12)and  (13)   developed  in  this  investigation.     Using  Ettle's  data,  these 
equations  produced  errors  in  equivalent  salinity  of    0.008°/oo.     Thus, 
equations  (12)  and  (13)  are  at  least  partially  verified,  but  more  eval- 
uation than  has  been  done  is  necessary. 

VI.      CONCLUSIONS  AND  RECOMMENDATIONS 

The  specific  conductivity  of  sea  water  at  temperatures  and 
salinities  representative  of  the  ocean,  increases  with  increasing 
pressure.    Mathematical  relationships  for  the  variations  indicate  that 
they  are  nonlinear.    The  most  common  method  of  previous  investiga- 
tions in  fitting  observed  data  has  been  by  using  polynomials  of  second 
or  third  degree.    However,  these  polynomials,  particularly  that  of 
Bradshaw  and  Schleicher ,  do  not  describe  the  relationship  of  con- 
ductivity and  pressure  adequately  at  pressures  near  atmospheric. 
Equations  (12  and  (13)  were  developed  relating  conductivity  ratio, 
temperature  and  pressure.    While  neither  equation  is  entirely  satis- 
factory ,  they  do  at  least  describe  the  present  data  and  the  data  of 
Ettle  well.    There  is  apparently  a  rather  complex  temperature-pres- 
sure effect  on  the  conductivity  of  sea  water ,  and  although  this  effect 
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may  be  small ,  it  can  be  significant ,  especially  if  accuracy  to  four 
significant  figures  is  required  for  in  situ  measurements  of  salinity. 
For  future  work,  more  detailed  comparison  of  equations  (12)  and 
(13)  developed  in  this  study  over  a  wider  range  of  salinities  is  recom- 
mended, and  more  data  at  temperatures  and  pressures  characteristic 
of  the  deep  ocean  would  be  useful.     Further  investigation  of  the  tem- 
perature-pressure effect  on  conductivity  is  also  recommended. 
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TABLE  XII.       Coefficients  of  the  Equation 

T-l  x  ln(R)  =  Ax  +  A2(P-PQ)  +  A3(P-P0): 
for  Sea  Water  of  Salinity  34.  099%o— 
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TABLE  XIII.     Coefficients  of  the  Equation 


(P-Po)"1  x  ln(R)= 

for  Sea  Water  of  Salinity  33.  643° /oo- 


B±  +B2T  +  B3T' 
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TABLE  XIV.     Coefficients  of  the  Equation 

(P-P  )"1x  ln(R)  =  Bi  +  B2T  +  B3T2 
for  Standard  Sea  Water  of 
Salinity  35.  001°/oo 
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TABLE  XV.       Coefficients  of  the  Equation 

(P-P0)-1  x  ln(R)  =  B1  +  B2T  +  B3T2 

for  Sea  Water  of  Salinity  34.  099°/oo 70 
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THIS  PROGRAM  TAKES  THE  RAW  EXPERIMENTAL  DATA  OF  TEMPERATURE, 
PRESSURE,  AND  PESISTANCE  AND  COMPUTES  THE  SPECIFIC  CON- 
DUCTANCE OF  SEA  WATER.   IT  WILL  ALSO  COMPUTE  THE  TEMPER- 
ATURE COEFFICIENTS  FOR  EACH  WATER  SAMPLE  AT  EACH 
PRESSURE  AND  THE  PRESSURE  COEFFICIENTS  AT  EACH 
TEMPERATURE. 

THE  VARIABLES  USED  IN  THIS  PROGRAM  ARE: 

SAL=SALINITY  OF  WATER  SAMPLES. 
P=PRESSURE  (PSI). 
PD=PRESSURE  (DECIBARS). 

RTN=RESISTANCE  OF  SAMPLE  "N1;  WHERE  ( =1 , 2 , 3 , . . . , K. 
CONDN=SPECIFIC  CONDUCTANCE  OF  SAMPLE  ■N'; 
WHERE  N=l,2,3,. . . ,K. 

THE  OTHER  VARIABLES  THAT  ARE  DIMENSIONED  ARE  USED  FOR 
TEMPORARY  STORAGE  OR  FOR  TITLING  PURPOSES. 

CONSTANTS  USED  IN  THIS  PROGRAM  ARE: 

K=NUMBER  OF  WATER  SAMPLES. 
L=NUMBER  OF  PRESSURES  MEASURED  (PSI). 
M=NUMBER  OF  TEMPERATURES  MEASURED  (DEG.  C). 
FACT=FACTOR  TO  CONVERT  PSI  TO  DECIBARS. 
DH=DIAMETER  OF  HYDRAULIC  PISTON  (INCHES). 
DV=DIAMETER  OF  PRESSURE  VESSEL  PISTON  (INCHES). 
CELL1=CELL  CONSTANT  FOR  CELL  ONE  (MHO/CM). 

CELL2=CELL  CONSTANT  FOR  CELL  TWO  (MHO/CM).   NOT  REQUIRED 
IF  ONLY  ONE  CELL  USED. 

REAL*8  ITITLEC 12,15), TITLE ( 10  ,  15 ) , DCON1 (  10)  ,DC0N2(10) , 
*DC0N3( 10),YAA(10,YAB(10),YAC(10),B(10),SB(10),DPD(10), 
*ITATLE( 12) ,TATLE(12),W(10) ,DELY( 10), 

REAL*8  ITOTLE( 1 2 , 10 ) , TOTLE ( 1 0  ,  10 )  ,C0N1( 15)  ,C0N2(15) , 
*C0N3( 15),YTA(15),YTB(15) ,YTC( 15) , C( 15 ) , SC ( 15 ) ,EPSY( 15) 
*,DT(15) ,WT( 15) 

REAL*8    DTRHK1  5)  ,DTRH2  (  15)  ,DTRH3(  15  )  ,  0XRH1  (  10  )  , 
*DXRH2( 10) ,DXRH3( 10) ,E( 15) ,DTRY(15) , XT A( 1 5 ) , SE ( 15 ) , 
*XTB(15,XTC( 15),D(10),SD(10),DERY(10),XAA(10),XAB(10), 
*XAC(10) 

REAL*8     ITETLE( 12 ),TETLE( 10) 

REAL*8    TTTLE(10)/10*'     •/ 

REAL  LABEL/4H     / 

DIMENSION  SAL (5) ,P( 10) ,PD(10),T(15) ,RT1( 15,10), 
*RT2(15, 10) ,RT3(15,10) ,C0ND1( 15,10) ,C0ND2( 15,10), 
*C0ND3(15,10) 

DIMENSION  DUM1 (10) ,DUM2( 10) ,DUM3( 10) ,TEM1(  10) ,TEM2( 10) 
*,TEM3( 10) 

DIMENSION  DMY1  (15) ,DMY2( 15) , DMY3( 15) ,TMP1(  15)  ,TMP2( 15) 
*,TMP3(15) 

DIMENSION  TRH1 (15) ,TRH2( 15) ,TRH3( 15) ,XRH1( 10) ,XRH2( 10) 
*XRH3(10) 

DATA  F A CT/O. 68 947/, OH/ 1. 875 /,DV/ 1.530/, CELL 1/0. 63 878/ 
*,CELL2/0. 75123/ 

READ  THE  RAW  DATA  IN. 

READ(5,1)  K,L,M 
1  F0RMAT(3I3) 
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READ(5,2)     (SAL ( I) ,1=1, K) 

2  F0RMAK3F6.3) 
READ(5,3)     (P(I),I=1,L) 

3  F0RMAT(9F5.0) 
READ(5,4)  (T(I),I=1,M) 

4  FORMAT (10F5. 2) 

READ(5,6)  (  (PTK  I,  J)  ,J=1,L)  ,I  =  1,M) 
RE AD (5,6)  ((RT2(I,J),J=1,L),I=1,M) 
READ(5,6)  (  (RT3(  I, J  ),J=1,L ),I=1,M) 

6  F0RMAT(9F6.3) 

READ  THE  INFORMATION  FOR  TITLING. 

READ(5,7)  (  (ITITLEC I, J)  ,I=1,12),J=1,M) 
READ(5,7)   ( (TITLE (  I,  J)  ,  1  =  1,10)  ,J  =  1,M) 

7  FORMAT (6A8) 

RE AD (5, 8)  ((IT0TLE(I,J),I=1,12),J=1,L) 
READ(5, 8)  (  (TOTLE( I  ,J)  ,1=1,10)  ,J  =  1,L) 

8  F0RMAT(6A8) 

ECHO  CHECK  THE  DATA. 

D010  1=1, K 

WRITE(6,1005)  SAL(I) 
1005  FORMAT( ' 1' , 5X,F6.3,/) 

D011  J=1,M 

WRITE(6,1015)  T(J) 
1015  FORMAT(5X,F5.2 ) 

IF(I.EQ.2)  GOTO  20 

IF(  I.EQ.3)  GOTO  30 

D012  N=1,L 

12  WRITE(6,1025)  P ( N ) , RT1 ( J ,N) 
GOTO  40 

20  D013  NN=1,L 

13  WRITE(6,1025)  P ( NN ) ,RT2 ( J ,NN ) 
GOTO  40 

30   D014  NM=1,L 

14  WRITE(6,1025)  P ( NM )  ,  RT3 ( J , NM ) 
1025  F0RMAT(5X,F5.0,5X,F6.3,/) 

40  CONTINUE 
11  CONTINUE 
10  CONTINUE 

CONVERT  PSI  TO  DECIBARS. 

D0100  I=1TL 
100  PD( I )=P( I J*( (DH/DV)**2)*FACT 

CONVERT  RESISTANCE  TO  OBSERVED  SPECIFIC  CONDUCTANCE. 

D0200  1=1, M 

D0210  J=1,L 

CONDK  I,  J)=1.0/RT1(  I, J) 

C0ND2( I, J)=1.0/RT2( I ,J) 

C0ND3( I, J)=1.0/RT3( I, J) 

REMOVE  THE  CARDS  FROM  HERE  THROUGH  GOTO  210  STATMENT  IF 
ONLY  ONE  CELL  HAS  BEEN  USED. 

IF(  I.GE.ll )  GOTO  250 

CONDK  I  ,  J)=C0ND1(  I,  J)*CELL2 

C0ND2( I , J)=C0ND2( I, J)*CELL2 

C0ND3(  I  ,  J)=C0ND3(I, J)*CELL2 

GOTO  210 
250  CONDK  I  ,  J)=C0ND1  (  I,  J)*CELL1 

C0ND2( I , J)=C0ND2( I, J)*CELL1 

C0ND3(  I  ,  J)=C0ND3( I , J)*CELL1 
210  CONTINUE 
200  CONTINUE 
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PRINT  THE  RESULTS. 

D0300  1=1, K 
WRITE(6,200P)  SAL(I) 
2000  F0RMAT( • 1* ,32X, '  SA LI NITY= • , 6X , F6. 3 , / ) 
D0310  J=1,M 
WRITE(6,2100)  T(J) 
2100  F0R,MAT(32Xt  '  TEMPERATURE=«  ,  3X  ,  F5  .  2  t  /  ) 

WRITP(6,2200) 
2200  F0RMAT(25X,«  PRESSURE ' , 2X , •  SPECIFIC  CONDUCTANCE',//) 
IFU.EQ.2)  GOTO  350 

IFU.EQ.3)  GOTO  360 
D0320  N=1,L 
320  WRITE(6,2300)  PD ( N ) , CONDI ( J , N ) 

GOTO  370 
350  D0330  NN=1,L 
330  WRITE(6,2300)  PD( NN ) , C0ND2 ( J , NN ) 

GOTO  370 
360  D0340  NM=1,L 

340  WRITE(6t2300)  PD ( NM  )  , C0ND3 ( J , NM ) 
2300  F0RMAT(26X,F6. 1,9X,F7.6,/) 
370  CONTINUE 
310  CONTINUE 
300  CONTINUE 

CONVERT  TO  DOUBLE  PRECISION  AND  TO  ONE  DIMENSION  ARRAYS  FOR 
SPECIFIC  CONDUCTANCE  AND  COMPUTE  THE  LEAST  SQUARES  POLY- 
NOMIAL FIT  TO  THE  DATA.   PLOT  THE  RESULTS  OF  THE  HIGHEST 
DEGREE  POLYNOMIAL  AND  THE  DATA.   (PRESSURE  COEFFICIENTS) 

D0400  1=1, M 

D0410  J  =  1,L 

DCONK  J  )=C0ND1  (  I»J) 

DC0N2( J )=C0ND2 ( I , J) 

DC0N3( J)=C0ND3 ( I , J) 

DPD( J)=PD( J ) 

DUMK  J)=C0MD1(  I,  J) 

DUM2( J)=C0N02( I, J) 

DUM3(J)=C0ND3( I, J) 

W( J)=l .ODO 
410  CONTINUE 

D0415  N=l,10 

TATLE(N)=TITLE(N,I ) 
415  CONTINUE 

THIS  SUBROUTINE  IN  NPGS  SOURCE  LIBRARY. 

CALL  LSQPL2(L,3,DPD,DC0Nl,W,YAA,DELY,BtSB,TATLE) 
CALL  LSQPL2(L,3,DPD,DC0N2,W,YAB,DELY,BtSB,TATLE) 
CALL  LSQPL2(L, 3 , DPD,DC0N3 T W , YAC , DELY, BT SB , TATLE ) 
D0420  NN=1,12 
ITATLE(NN)=ITITLE(NN,I ) 
420  CONTINUE 

THE  0.025  IS  A  SCALING  FACTOR  FOR  THE  DATA  OF  THIS  STUDY. 

D0425  NM=1»L 

TEM1(NM)=YAA(NM)-0.025D0 
TEM2(NM)=YAB(NM)-0.02  5D0 
TEM3(NM)=YAC(NM)-0.025DO 
DUM1(NM)=0UM1( NM) -0.02  5 
DUM2(NM )=DUM2(NM)-0.02  5 
DUM3(NM)=DUM3(NM)-0.02  5 
425  CONTINUE 

THIS  SUBROUTINE  IN  NPGS  SOURCE  LIBRARY. 

CALL  DRAW(L,PD tDUMl t 1,1, LABEL, ITATLE , 70  0. 0 , 0. 003 , 0, 0, 0 
*, 0,6, 8,0, LAST) 

CALL  DRAW(L,PD ,TEM1 ,3 , 0 , LABEL , ITATL E, 700 . 0, 0 .003 , 0, 0, 0 
*, 0,6, 8,0, LAST) 
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CALL  DPAW(L,PD,DUM2,1,2,LABEL, I  TATL E, 70  0 . 0,0  .  003 , 0, 0,  0 
*, 0,6, 8,0, LAST) 

CALL  DR AW (L,PD,TEM2, 3,0, LABEL, IT ATL E , 700 . 0 , 0 .003 , 0, 0, 0 
*,0,6,8,0,LAST) 

CALL  DRAW(L,PD,DUM3,1,3,LABEL, IT ATL E , 700 . 0 , 0 . 003 , 0, 0 , 0 
*,0,6,8,0,LAST) 

CALL  DRAW(L,PD,TEM3,3,0,LABEL,  IT ATLE , 700 . 0, 0 . 003,  0,  0,  0 
*, 0,6, 8,0, LAST) 

400  CONTINUE 

CONVERT  TO  DOUBLE  PRECISION  AND  TO  ONE  DIMENSION  ARRAY  FOR 

SPECIFIC  CONDUCTANCE  AND  COMPUTE  THE  LEAST  SQUARES  POLY- 
NOMIAL FIT  TO  THE  DATA.  PLOT  THE  RESULTS  OF  THE  HIGHEST 
DEGREE  POLYNOMIAL  AND  THE  DATA.   (TEMPERATURE  COEFF.) 

D0401  J=1,L 

D0411  I=1,M 

C0N1 ( I )  =  C0ND1(  I,  J) 

C0N2( I )=C0ND2(  I, J) 

C0N3( I )=C0ND3(  I,  J) 

DT( I )=T(I) 

DMYK  I  )=C0MD1(  I, J) -0.025 

DMY2( I )=C0ND2( I, J) -0.02 5 

DMY3U  )=C0ND3(  I,  JJ-0.025 

WT( I )=1 .ODO 

411  CONTINUE 
D0412  N=l,10 
TETLE(N)=TOTLE (N,J) 

412  CONTINUE 
KM=3 

KM=THE  DEGREE  OF  FIT. 

CALL  LSQPL2(M, KM, DT ,C0N1 , WT, YT A, EPSY, C, SC , TETLE ) 
CALL  LS0PL2(M, KM, DT,CCN2,WT,YTB, EPSY, C,SC, TETLE) 
CALL  LSQPL2(M, KM , DT ,C0N3 , WT ,YTC, EPSY, C , SC , TETLE ) 
D0413  NN=1,12 
ITETLE(NN)=ITOTLE(NN,J) 

413  CONTINUE 
D0414  NM=1,M 

TMPKNM  )=YTA(NM)-0.025D0 
TMP2(NM)=YTB(NM) -0,025 DO 
TMP3(NM)=YTC(NM)-0.025D0 

414  CONTINUE 

CALL  DRAW(M,T,DMY1,1  ,1 ,LABELt I TETLE, 4.0 , 0 . 003 ,0, 0 ,0 , 0, 
*6, 8,0, LAST) 

CALL  DRAW(M,T,TMP1 ,3,0, LAB EL, I TETLE , 4.0 , 0 . 003 ,0, 0 , 0, 0, 
*6, 8,0, LAST) 

CALL  DRAW(M,T,DMY2,1,2 , LAB  EL, I TETLE ,4.0 , 0. 003 , 0 , 0, 0, 0, 
*6,  8,0, LAST) 

CALL  DR AW (M,T,TMP2, 3,0, LABEL, I TETLE , 4.0 , 0. 003 , 0,0 , 0, 0, 
*6, 8,0, LAST) 

CALL  DRAW(M,T,DMY3, 1,3, LABEL, I TETLE , 4.0 , 0 . 003 , 0 ,0,0, 0, 
*6, 8,0, LAST) 

CALL  DRAW(M,T,TMP3,3 ,0, LAB  EL, I TETLE ,4.0 , 0 . 003 ,0, 0, 0, 0, 
*6, 8,0, LAST) 

401  CONTINUE 

COMPUTE  THE  LEAST  SQUARES  FIT  TO  THE  OBSERVED  CONDUCTIVITY 
RATIO.   THE  FORM  OF  THE  EQUATION  DESIRED  IS: 
LN(R)/(P-P1)=A(1 )+A(2 )*T  +  A(3  )*T**2 

D0167    1=1, M 

T(I  )=T(  D+273.16 

DT( I)=T(  I) 
167  CONTINUE 

PD(1)=10.13 

D0137  J=1,L 

PD( J)=PD( J)-PD(l) 

DPD( J)=PD(J) 
137  CONTINUE 
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168 


169 


WT( 1)=1.0D0 
D0169    J=2tL 
D0168    I =1 , M 
TRH1 (I )=ALGG(RH1( I 
TRH2(I )=AL0G(RH2( I 
TRH3< I )=AL0G(RH3( I 
WT( I )=1.0D0 
DTRHK  I  )=TRH1(  I) 
DTRH2( I )=TRH2( I ) 
DTRH3( I )=TRH3( I ) 
CONTINUE 
DTRHK  1  )=0.0D0 
DTRH2( 1 )=0.0D0 
DTRH3( 1 )=0.0D0 
CALL    LSQPL2(M,3 
CALL    LSQPL2(M,3 
CALL    LSQPL2(M,3 
CONTINUE 


DT, 
DT, 
DT, 


) )/PD( J 
) )/PD( J 
) )/PD(J 


DTRHltWT,XTATDTRY,E,SE,TTTLE) 
DTRH2,WT,XTB,DTRY,E,SE,TTTLE) 
DTRH3,WT,XTC,DTRY,E,SE»TTTLE) 


COMPUTE  THE  LEAST  SQUARES  FIT  TO  THE  OBSERVED  CONDUCTIVITY 
RATIO.   THE  FORM  OF  THE  EQUATION  DESIRED  IS: 
LN(R)/T=A(1 )+A(2 )*(P-P1 )+A(3 )*(P-P1)**2 


198 


199 


D0199  I=1,M 

D0198  J=1,L 

XRHK  J)=AL0G(RH1(  I 

XRH2( J)=AL0G(RH2( I 

XRH3(J)=AL0G(RH3( I 

DXRHK  J)=XRH1(  J) 

DXRH2( J)=XRH2( J) 

DXRH3( J)=XRH3( J) 

W( J)=1.0D0 

CONTINUE 

CALL  LSQPL2(L,3 

CALL  LSQPL2(L,3 

CALL  LSQPL2(L,3 

CONTINUE 

STOP 

END 


J)  )/T(I ) 
J) )/T(  I  ) 
)/T( 


J) 


I  ) 


DPD,DXRH1TW,XAA,DERY,D,SD,TTTLE) 
DPD,DXRH2f W,XABtDERY,D,SDf TTTLE ) 
DPD,DXRH3,W,XAC,DERY,D,SD,TTTLE) 
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